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We develop an approach to use nanostructured plasmonic materials as a 
non-magnetic negative-refractive index system at optical and near-infrared 
frequencies. In contrast to conventional negative refraction materials, our 
design does not require periodicity and thus is highly tolerant to fabrication 
defects. Moreover, since the proposed materials are intrinsically non-magnetic, 
their performance is not limited to proximity of a resonance so that the re- 
sulting structure has relatively low loss. We develop the analytical description 
of the relevant electromagnetic phenomena and justify our analytic results via 
numerical solutions of Maxwell equations. © 2008 Optical Society of America 

OCIS codes: 999.9999 Negative index of refraction, 160.4760 Optical properties, 
110.2990 Image formation theory, 350.5730 Resolution 

1. Introduction 

The electromagnetic response of negative refractive index materials (NIM^P^^EEEEIEl h as 
recently attracted unprecedented attention. Novel optical phenomena, predicted to take place 
in these unique systems include reversal of Snell Law, Doppler Effect, Cherenkov Effect P 
aberration-free^^ and sub-difrractiorl ^ * 1 ^ 11 * 12 * 1 ^ 1 ^ 1 ^ 1 ^ imaging, and excitation of the new 
types of surface and nonlinear waves.^^ In particular, realization of NIMs may potentially 
lead to fabrication of new types of lenses and prisms^^ new lithographic techniques,' 9 -'^' 
novel radars, sensors, and telecommunication systems. However, despite the great advantages 
NIM has to offer for optical and infrared spectral range, all practical realizations of NIM are 
currently limited to GHz frequencies.' 19 * 20 * 2 *^' 
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Until recently there were two major approaches for NIM design. The first one is based 
on the original proposaP that material with simultaneously negative dielectric permittivity 
and magnetic permeability must have a negative refraction index. This particular approach 
also benefits from the possibility to resonantly excite the plasmon polariton waves at the 
interface between NIM and surrounding media, which in turn may lead to sub-diffraction 
imaging) 3 * 9 * 10 * 11 * 12 * 13 * 14 * However, the absence of natural magnetism at high (optical or in- 
frared) frequencieipSl requires the design and fabrication of nanostructured meta-materials, to 
achieve the non-trivial magnetic permeability.*^ ^ 4 * 25 * 26 * 27 * 28 * 29 * As these engineered systems 
typically operate in close proximity to resonance, resonant losses become the dominant factor 
in system response, severely limiting the practicality of resonant-based systems! 12 * 13 * 11 * 10 ^* 

The second approach for NIM design involves the use of photonic crystals! 21 * 32 * 33 * 22 -*^ How- 
ever, the NIM response in these systems is typically associated with second or other higher- 
order bands and requires a complete bandgap between the band in use and the next band. 
The dispersion and very existence of the required bandgap are typically strongly affected by 
crystal disorder, unavoidable during the fabrication step. The manufacturing of the optical 
photonic crystals-based NIM typically requires 3D patterning with 10-ran - accuracy, which 
is beyond the capabilities of modern technology. 

To address the mentioned- above shortcomings of the traditional NIM schemes, we have 
recently introduced an alternative approach to design the NIM structure.^ In contrast to 
"conventional" systems, the proposed design does not rely on either magnetism or periodicity 
to achieve negative refraction response. It has been shown that the combination of strong 
anisotropy of the dielectric constant and planar waveguide geometry yields the required 
negative phase velocity in the system.^ Here we present the detailed description of NIMs 
proposed in Refs. 0011, study the effects related to waveguide boundaries, important for 
optical domain, and suggest several nanostructured materials providing the low-loss negative 
refraction response at optical and infrared frequencies. 

The rest of the paper is organized as follows: the next Section is devoted to EM wave 
propagation in strongly anisotropic waveguides; Section 3 describes the proposed realizations 
of the structure; imagining properties of these composites are shown in Section 4; Section 5 
concludes the paper. 

2. Negative refraction in strongly anisotropic waveguides 

We consider wave propagation in the 2D planar waveguide structure shown in Fig. ^ The 
propagation in the system is allowed in the y and z directions, while the waveguide walls 
occupy the regions \x\ > d/2. The waveguide core is assumed to be a homogeneous, non- 
magnetic (fi = 1) material, with a uniaxial anisotropic dielectric constant with dielectric 
permittivities e± and en along and perpendicular to the optical axis respectively. The optical 
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axis of the core material (C) is assumed to be perpendicular to the direction of the wave 
propagation in the media (C||x). Therefore, despite the anisotropy of the system, the effective 
refractive index of propagating in the planar geometry waves will be completely isotropic. 




Fig. 1. Schematic configuration of non- magnetic negative-refraction system 



Any wave propagating in the system can be represented as a linear combination of the 
waveguide modesW^ An individual mode is defined by its structure along the optical axis 
direction (C) and its polarization. Two different kinds of modes have to be distinguished. 
The modes of the first kind (known as TE waves) have their E vector perpendicular to 
the optical axis. The propagation of these waves is fully described by the in-plane dielectric 
constant ey. The modes of the second kind (known as TM waves) have their H vector in 
the waveguide plane and are affected by both e±_ and ey. The existence of these TM waves is 
crucial for the NIM described here. 

In the analytical results presented below we limit ourselves to the case of single-mode 
propagation. We note that such a description provides complete information about the linear 
properties of the waveguide structure. Indeed, as mentioned above, an arbitrary wavepacket 
in the system can be represented as a linear combination of modes. In our numerical simu- 
lations discussed in Section 4 we utilize this property to compute the imaging performance 
of the system. 

2. A. Waveguide with perfectly conducting walls 

As it has been shown in Ref. jjj, the propagation of a mode in a planar waveguide can be 
described by the free-space-like dispersion relation: 

k 2 y + k] = eidc 2 , (1) 

where e is en for TE modes and e± for TM ones, k y and k z are the propagation components of 
the wavevector, k = uj/c with a> and c being the free-space angular frequency of the radiation, 
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and speed of light in a vacuum; the propagation constant v is given by 
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and the parameter k defines the mode structure in x direction. 

As it directly follows from Eq. (JTJ), the phase velocity of a propagating mode is equal to 

v p = n k, (3) 

where the effective refraction index n 2 = ev. Note that similar to the case of the plane wave 
propagation in free-space, the refraction index contains a product of two (mode-specific) 
scalar constants. A transparent structure must have both propagation constants of the same 
sign. The case of positive e and v corresponds to "conventional" (positive refraction index) 
material. The case of negative e and v describes NIM 12124] The NIM behavior can be easily 
illustrated by comparing the Poynting vector S z and the wavevector k z as shown below. 

Similar to any waveguide structure, the mode in the system described here can be re- 
lated to the x profile of the longitudinal field component (the detailed description of such a 
dependence is given in Ref. [Zj). To better illustrate the physical picture behind the mode 
propagation, in this section we present the results for the important case of perfectly conduct- 
ing waveguide walls. In this case, the EM energy is confined to the waveguide core and the 
longitudinal field has a cos(nx) or sin(nx) profile depending on the symmetry with respect 
to the x = plane, with k = (2j + l)7r/d for symmetric and k = 2irj/d for anti-symmetric 
modes respectively with j being the integer mode number. The deviation from this idealized 
picture due to finite conductance of the waveguide material does not change the physical 
picture described in this section, and for the practical case of "good" metals (Ag,Al,Au) at 
near-IR to THz frequencies can be treated perturbatively. Results of such a perturbation 
approach are presented in the Section 2b. 

The electric (Ue) and magnetic (Uh) field contribution to the energy density of a mode 
in weakly-dispersive material (\e/u\ |de/da>|) is given by Ue = j^b /(D • E*)dx and 
Uh = J(H ■ H*)dx respectively^ (the asterisk (*) denotes the complex conjugation). 
Using the explicit mode structure for TE and TM waves (see Ref. [2j) we arrive to: 

2l 2 2l 2 

jj(jTM) = jj(TM) = _l_fit| Ao | 2; U(TM } = U{ TM) + ^M, = ^j^. (4) 

ur ] = ur ] = ^oi 2 ; = ur + ur ] = ^ 

WIT 07T 

where Aq is the mode amplitude. Thus, extending the similarity between the waveguide 
system described here and the free-space propagation, the EM energy of any propagating 
wave is always positive and contains equal contributions from the electric and magnetic 
components of the field. 
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It is also seen that the TE mode is in some sense very similar to the conventional plane 
wave propagating in the isotropic homogeneous dielectric. Namely, (i) energy density of the 
TE waves is exactly equal to that of the plane waves; (ii) there is no wave propagation in 
material with en < 0. In contrast to this behavior, the sign of the dielectric permittivity 
alone does not impose limitations on the propagation of TM modes. 

Another important characteristic of the energy transport in the EM system is the average 
energy flux given by the propagating component of the Poynting vector S = ^-[E x H]. 
Selecting the direction of the wave propagation as z axis, we obtain: 

S ({TE; TM}) _ c k * jj(\TE; TM}) ^ 

e{||;±} k 

It is clearly seen from Eq. El that the relation between the direction of the phase velocity 
and energy flux is defined by the sign of the dielectric constant (for a given mode polarization) 
- positive e means n > propagation, while e < signifies the NIM case. Of course, for this 
relation to take place, we must require the medium to be transparent - both propagation 
constants e and v must be of the same sign. As it can be seen from Eq. (Q), the NIM condition 
can be satisfied only for TM wave and only in the case of extreme anisotropy of the dielectric 
constant of the core material (eii ej_ < 0). The feasibility of the fabrication of such unusual 
materials will be discussed in the Section 3. 

2.B. The effect of finite wall conductance 

In this section we consider the practical realization of the system described above, in which 
the anisotropic core material is surrounded by metallic walls. The electromagnetic properties 
of metals at high (GHz to optical) frequencies are dominated by the dynamics of the free- 
electron plasma-like gas. Following the approach described in e.g. Ref. [33] it is possible to 
write down the high-frequency effective permittivity of metal in Drude form: 

e m (u) = — — — , (7) 

U)(u + ir) 

where the constant term eoo describes the contribution of the bound electrons, r is responsible 
for EM losses due to (inelastic) processes, and Q p \ = is the plasma frequency with N e , e, 
and ra e ff being the free-electron concentration, charge, and effective mass respectively. Note 
that for to < f2 p i/\/^°° the permittivity of the metal becomes negative e' m < (here and 
below single and double prime ('; ") denote the real and imaginary parts respectively). For 
most of "good" metals (Ag,Al,Au) the plasma frequency is of the order of 10 eV and ~ 1, 
which means that e' m is negative from optical to GHz frequencies. The losses, given by the 
parameter e^/|eJJ <C 1 are typically small in these spectral ranges. 

Similar to the case of perfectly conducting waveguide walls, the structure of the modes in 
the system can be still derived from the dependence of the longitudinal (z) field component 
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on the x coordinate, which has cos(kx) or sin(Kx) behavior depending on its symmetry. The 
exact value of the mode parameter k is given by the requirement of the in-plane (y, z) field 
components continuity throughout x = ±d/2 planes. For symmetric (cos) mode profile, we 
obtain: 




In the limit of e m — > — oo, these equations yield the values k,q = ir(2j + l)/d, used in 
the previous Section. As we previously noted, these values correspond to the well-known 
condition of zero mode magnitude at the waveguide boundary. In the limit of sufficiently 
large \e m \ it is possible to find the correction to the above values of the mode parameter k. 
Specifically, 

k<™ ~ K Jl-^L=\ (9) 

K {TE) _ Ko(l- 2 ^ 

As the mode parameter k plays a role of an inverse confinement length of the mode in x 
direction, the negative k correction signifies the "mode expansion" into the waveguide wall 
region. Such a mode expansion is illustrated in Fig. |2| 

E„AU E ( ,AU 




Fig. 2. The cross-section of the field in the planar waveguide with hollow d = 
0.5 /im-thick core; dashed line corresponds to the case of e m = — oo (perfect 
metal boundary); solid line corresponds to Ag boundaries for A = 0.85 fim; 
dots correspond to k calculated using Eq. ©; TM (a) and TE (b) modes are 
shown 



The immediate effect of such a change in the mode structure is the change of the effective 
phase velocity, given by the refraction index: 

n<™> « ( 1 + w 2 / ) (10) 
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where vq = 1 — KQ/(e||k 2 ). As it has been described above, the sign of the refraction index 
for the TM polarization has to be selected positive for e± > 0; v > 0, and negative for 
e_L < 0; v < 0. 

Penetration of the mode into the waveguide wall region has another effect on the wave 
propagation. Namely, the finite value of the e" m introduces an additional (with respect to 
the core material) absorption into the system. As a result, the magnitude of a mode will 
exponentially decay as it propagates through the system. Such an attenuation can be related 
to the imaginary part of the effective refractive index through E oc exp(— n"kz). In the limit 
of small absorption in the metal (e^/|e^J <C 1) the "waveguide-induced" mode decay is 
described by: 



regardless of the sign of the refractive index. 

Using Eq. (jllj) we estimate that for wavelengths A > 850 nm, the losses introduced by a 
silver waveguide walls are substantially small (n"/n < 0.01). 

3. Anisotropic nanoplasmonic composites 

We now consider the fabrication perspectives of the material with strong optical anisotropy 
required for NIM waveguide core region. A number of naturally occurring materials with 
the required properties exist at THz or far IR frequencies. Some examples include Bi and 
Sapphire.^ Unfortunately, no known material exhibits anisotropy exceeding 30% at optical or 
infrared spectral range. Here we propose to take advantage of a new class of nano-engineered 
media, known as meta-materialsP^ In these composites, nanostructured particles are used as 
meta-atoms to achieve the desired EM properties. 

To realize the strong optical anisotropy we propose to use a combination of plasmonic or 
polar particles (providing the negative permittivity) and dielectric media (having e > 0). If 
the characteristic size of inhomogeneities and their typical separation are much smaller than 
the wavelength of incident radiation, the EM response of the composite structure can be 
described in terms of the effective dielectric constant e e gW^ 



k 3 d^e\\u \e m \ \e' m \ 

Note that in agreement with causality principle^* 2 ^ the losses in the system are positive, 




< D(r) > a =< e{r) a>p E{r)p >= e eSaf! < E(r) > p , 



(12) 
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where the brackets (<>) denote the averaging over the microscopically-large (multi-particle), 
macroscopically small (subwavelength) spatial area, Greek indices denote Cartesian compo- 
nents, and summation over repeated indices is assumed. Since the size of a particle a typically 
enters Maxwell equations in the combination kaj^ all size-related effects play a minor role 
in the considered here "quasi-static" averaging process. Therefore, we note that the compos- 
ites proposed below are highly tolerant with respect to size variation. Also, since the effects 
described here are obtained in effective medium approximation, the desired response does 
not require any periodicity of the particle arrangement and only the average concentration 
has to be controlled during the fabrication step. 

Below we present two meta-material designs of the strongly anisotropic composite for 
optical and infrared spectrum ranges. 

3. A. Layered system 

We first consider the permittivity of a stack of interlacing plasmonic (Ag,Au,Al,. . .) or po- 
lar (SiC) (e p i < 0) and dielectric (Si,GaAs,. . .) (e^ > 0) layers. We assume that the layers 
are aligned in the waveguide (y, z) plane (see Fig. EJ). As noted above, the absolute thick- 
ness of the layers is not important (as long as it is subwavelength), and only the average 
concentration of plasmonic layers N p i plays the role. 




Fig. 3. The schematics of the layered structure described in the text 



To compute e e fr we note that the E y ,E z , and eE x have to be continuous throughout the 
system,^ 1 41 ' 42 1 leading to: 

£|| = e e fr B , z = N p i e pi + (1 - N pl )e d (13) 

6± ~ €eSx ~ (1-N pl )e pl + N pl e d 

The effective permittivities for several layered composites are shown in Fig. 0] We note 
that while the strong anisotropy ey • < can be easily achieved in the layered system, the 
actual realizations of the materials with en > 0, e± < required for the high-frequency NIM 
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Fig. 4. Real part (a,c,e) and absorption (b,d,f) of effective e± (solid lines) and 
e|| (dashed lines) for layered systems; (a,b): Ag — Si stack; N p \ = 0.6; (c,d): 
Ag — SiOi stack; N p \ = 0.1; Note the extremely small absorption of this system; 
(e,f) SiC-Si stack; iV pl = 0.1. 



described here typically have substantial absorption 1 , and therefore have a limited range of 
applications.^^ 

On the contrary, the materials with en < 0, > (achieved, for example by a repeated 
deposition of Ag-Si layers) form low-loss media. While this configuration has a positive refrac- 
tion index, it may be potentially used to concentrate propagating modes in subwavelength 
areas.^^ 

3.B. Aligned wire structure 

The array of aligned e p i < nanowires embedded in the dielectric (e^ > 0) host, schematically 
shown in Fig. is in some sense a counterpart of the layered system described above. In fact, 
the boundary conditions now require the continuity of the E x field, along with the solution of 
quasi-static equations in the y — z plane. While in the general case the analytical solution of 
this problem is complicated, the case of small plasmonic material concentration is adequately 
described by the Maxwell-Garnett approximat ion p B l ?7 l ?$ l ?9 l 

_ N p i e p iE in + (1 — N p i)edE 

6 " ~ eeffy " ~ N pl E in + (1 - N pl )E [ ' 

e± = fieff, = N p i€ p i + (1 - N p i)e d , 

1 This particular realization of layered NIM structure for IR frequencies has been earlier proposed in 
Rcf. g] 
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where E in = _*; d E is the field inside the plasmonic inclusion and E is the excitation field. 




a) b) c) 



Fig. 5. (a) The schematics of the wired structure described in the text, (b- 
c) comparison of < e >||, calculated using Eq. (f*H|) (solid line) and derived 
from numerical solution of Maxwell Eqs as described in the text (dots); de- 
pendence of < e >|| on dielectric constant of the inclusions for N p \ (b) and on 
concentration for e p i = — 10 (c) is shown. 

To illustrate the validity of the MG approximation, we numerically solve the Maxwell 
equations in the planar geometry using the coupled- dip ole approach (CDA), described in 
detail in Refs. )2*%ll39ll40| . In these calculations the composite is represented by a large number 
of interacting point dipoles, and the resulting dipole moment distribution is related to the 
effective dielectric constant. Fig. 03 shows the excellent agreement between the numerical 
simulations and the analytical result [Eq. 1)14)1] . 

The effective dielectric constants for some composite materials are presented in Fig. |H1 
Note that in contrast to the layered system described above, these wired composites have 
extremely low absorption in the near-IR-NIM regime - in a way solving the major problem 
with the "conventional" design of optical LHMs. 

4. Imaging Properties of non-magnetic optical NIMs 

To illustrate the imaging performance of the proposed system we calculate the propagation 
of a wavepacket formed by a double-slit source through the 5//m-long planar layer of 5% 
Ag, 95% Si0 2 wire-based NIM-core described above [see Fig. Efa,b)], embedded in the Si 
waveguide. We select the thickness of the dielectric core to be d = 0.3/im, and assume the 
excitation by the telecom-wavelength A = 1.5/im. The Eqs. f)10pi4jl yield the following values 
of the refraction index: « 2.6, rS LHM ^ ps —2.6 + 0.05i 

To calculate the resulting field distribution we first represent the wavepacket at the z = 
plane as a linear combination of the waveguide modes.^^ We then use the boundary 
conditions at the front and back interfaces of the NIM region to calculate the reflection and 
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Fig. 6. Real part (a,c) and absorption (b,d) of effective e± (solid lines) and e\\ 
(dashed lines) for wired systems; (a,b): Ag — SiC>2 structure (note the relatively 
small absorption for the NIM regime); N p \ = 0.05; (c,d): SiC — Si structure; 
7V pl = 0.1. 



transmission of individual mode. The solutions of Maxwell equations are then represented 
as a sum of solutions for the individual modes. 

To better illustrate the imaging properties of the system and distinguish between the effects 
of negative refractive index and material absorption, we first neglect losses in the NIM core. 
The resulting intensity distribution in the system is shown in Fig.Efa). The image formation 
in the focal plane (z = 10/im) of the far-field planar NIM lens is clearly seen. 

In Fig. Wtjo) we compare the imaging through the planar NIM lens with and without the 
material absorption and demonstrate that the presence of weak loss, although it reduces the 
magnitude of the signal, it does not destroy the far-field imaging. Similar to any far-field 
imaging system ffi 1 1 ^^ 1 the resolution A of the non-magnetic NIM structure presented here 
is limited by the internal wavelength: A « Aj„/2 = A/(2n) ~ 0.3/im [see Fig. E^c)]. 

5. Conclusions 

In conclusion, we presented a non-magnetic non-periodic design of a system with negative in- 
dex of refraction. We have further proposed several low-loss nanoplasmonic-based realizations 
of the proposed structure for optical and infrared frequencies. We have presented analyti- 
cal description of the effective dielectric permittivity of strongly anisotropic nanostructured 
composites, and showed the excellent agreement of the developed theory with results of nu- 
merical solution of Maxwell equations. Finally, we have demonstrated the low-loss far-field 
planar NIM lens for A = 1.5/xm with resolution A « 0.3/im. 

The authors would like to thank E.E. Mishchenko and A.L. Efros for fruitful discussions. 
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Fig. 7. Imaging by a planar NIM-based lens, n > region: S^-filled planar 
waveguide; d = 0.3/xm; NIM region: planar waveguide with core material de- 
scribed in Fig. Efa-b); (a) The intensity distribution in the system with ab- 
sorption losses neglected; LHM region is beween z = 2.5/xm and z = 7.5/xm; 
focal plane corresponds to z — 10/im (white dashed line); slit size w = 0.75/im; 
(b) dashed line: emitted radiation; solid line: focal plane intensity distribution 
in system described in (a); dash-dotted line: same as solid line, but in the case 
of real (absorbing) NIM. (c) same as (b), but w = 0.3/xm (corresponding to 
far-field resolution limit of the system) 



The research was partially supported by NSF grants DMR-0134736, ECS-0400615 and Ore- 
gon State University. 
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